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We show that a nonlinear Hamiltoiiian evolution can transform an SU(3) coherent state into a 
superposition of distinct SU(3) coherent states, with a superposition of two SU(2) coherent states 
presented as a special case. A phase space representation is depicted by projecting the multi- 
dimensional Q-symbol for the state to a spherical subdomain of the coset space. We discuss realiza- 
tions of this nonlinear evolution in the contexts of nonlinear optics and Bose-Einstein condensates. 



INTRODUCTION 
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The superposition principle in quantum physics implies that superpositions of probability amplitudes for classical- 
like states are possible, yet superpositions of quasi-classical quantum field states are not generally observed in practice. 
Such superposition states are especially interesting as they dramatically illustrate the quantum superposition principle 
by creating coherent superpositions of distinct states, with these distinct states behaving as classical physical states. 
Extensive theoretical and experimental studies have been directed towards understanding and obtaining superpositions 
of distinct quasi-classical states, often referred to as Schrodinger cat states jl], |). In general the interest has been in 
generating superpositions of Heisenberg-Weyl coherent states via a nonlinear Hamiltonian evolution Jl], || , where the 
Heisenberg-Weyl coherent states are the displaced harmonic-oscillator vacuum state, 
j-^. Heisenberg-Weyl coherent states are relevant for studying the quantum-classical transition for the harmonic oscilla- 

• tor with one degree of freedom. The generalization to harmonic oscillators with two or more degrees of freedom leads 
to superpositions of product coherent states, or entangled coherent states, with a wealth of phenomena to be studied 
for such systems[Q. However, another generalization of superpositions of coherent states is possible by employing the 
^ generalized coherent states for general group actionsjH). Examples include superpositions of SU(2) coherent states|| 
■ and superpositions of SU(1,1) coherent statesQ, as well as superpositions of multiparticle SU(2) and SU(1,1) coherent 
states^. Whereas the Heisenberg-Weyl group is the symmetry group for harmonic oscillators, SU(2) is the symmetry 
group for spin precession, for two-channel interferometry and for the dynamics of ideal two-level atoms[^] ^] [To] pd| , 
and SU(1,1) is the symmetry group for the production of squeezed light and quantum interferometry with parametric 
up- and down-conversion |t], pT[ . The generalized coherent states for such systems represent the quasi-classical states. 
For example the SU(2) coherent states are also known as atomic coherent states and are analogous to classical elec- 
^ | trie dipoles|]lO|. Superpositions of these generalized coherent states in quantum systems are therefore of interest in 
studying counter-intuitive manifestations of the superposition principle. 

The studies of superpositions of Heisenberg-Weyl, SU(2) and SU(1,1) coherent states are simplified by the fact that 
these are one-parameter groups, and the corresponding coset space is thus a locally two-dimensional manifold. For 
the Heisenberg-Weyl group, the manifold is the plane, for SU(2) the Poincare sphere and for SU(1,1) the Lobachevsky 
plane. Here we generalize the studies of superpositions of coherent states to a two-parameter Lie group, namely SU(3). 
The manifold for phase space dynamics is the four-dimensional space isomorphic to the coset space SU(3)/U(2). 
Working in a manifold which is greater than two dimensions presents problems in terms of visualizing the dynamics, 
which we resolve by projecting the dynamics to subdomains of spherical manifolds. 

In order to generate superpositions of SU(3) coherent states, we consider the three-dimensional nonlinear oscillator. 
Such a model can be realized by employing the three-boson realization of the SU(3) generators and constructing the 
dynamics optically through employing passive optical devices and Kerr-type nonlincarities. Another physical system 
involves three interacting, independent Bose-Einstein condensates (BECs) []l2]. In order to bring out the essential 
properties of superpositions of SU(3) coherent states, we treat the dynamics as a closed system. 

Superpositions of SU(3) coherent states may be interesting from the perspective of weak force detection and accurate 
measurements. For example, weak force detection for two coupled, separated BECs employs a superposition of two 
extremal SU(2) coherent statcsjl3]. Superpositions of two extremal SU(2) coherent states are valuable for related 
applications, such as precision measurement in particle interferometry jl4|. Symmetries of SU(3) coherent states may 
also prove to be valuable in the context of weak force detection for three coupled BECs or for three-channel particle 
interferometry, as examples. 

We show that an initial SU(3) coherent state in the nonlinear oscillator evolves to a superposition of SU(3) coherent 
states. An analytical solution is provided which shows explicitly that a superposition of sixteen SU(3) coherent states 
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occurs. The analytical solutions for the SU(3) case are more complex and interesting than those for SU(2), and we 
elaborate on the superposition of two SU(2) coherent states as a special case. Graphical solutions are employed by 
using the Q-symbol (sometimes referred to as the Q-function) for the state and projecting to spherical subdomains. 
The graphical method provides an intuitive picture of the dynamics. 

THREE-DIMENSIONAL NONLINEAR OSCILLATOR 

The Hamiltonian for the three-dimensional nonlinear oscillator includes three pairs of annihilation and creation 
operators, designated by c\ , c\ , £2 , c| , £3 , £3 , with each pair {ci,c\} corresponding to one degree of freedom, or "mode" . 
The usual boson commutation relation [cfe,ct] = Skj applies. We also define number operators as = c\,Ck, and a 

total number operator N — J2 k=1 rife. We consider an isolated system; hence, total energy and total particle number 
are conserved. In particular we are interested in the simplest system with energy and number conservation which will 
allow an SU(3) coherent state to evolve into a superposition of SU(3) coherent states. Such a Hamiltonian is given 
by 

33 33 
H = uj^h k + xi(t) y^»j(^fc - 1) + Xs(t) ^2 fikfij + ^2 ^jk(t)CjC k , (1) 

fe=l fe=l j,k = l j,k=l 

j=£k j^k 

with Qjk(t) — Q k j(t). In this model, the coupling constant xi corresponds to the strength of the self-modulation 
term, and %2 is the strength of the cross-modulation term. The time-dependent quantities Xi are assumed to be 
controllable as well as the time-dependent coupling strengths Qjf. that are responsible for linear coupling between 
modes. For the case that each mode corresponds to bosons in a particular spatial region, the terms Xi quantify the 
nonlinear interactions, and Cljk quantify the strength of (linear) quantum tunneling of bosons between modes. For 
Xi = 0, with i = 1,2, the Hamiltonian (|lj) is a linear combination of su(3) generators, and thus the evolution is linear. 

The Hamiltonian (|l|) corresponds to two potential physical realizations, the first being a nonlinear optical four- 
wave mixing medium, that is a medium with a x^ nonlinearity. In four-wave mixing, the three boson operators c k 
correspond to single-mode fields interacting in the medium, and the coefficients %i, * € {1,2} are obtained by 
judiciously constructing the x tensor coefficients. 

A second realization arises in the context of interacting BECs with a small number of particles Jl2[ which is of 
interest for some condensates fl6|| . The three boson operators correspond to localized BEC modes with some overlap 
between the modes (in contrast to the case considered in Ref. [[l2| , which assumes a negligible overlap) . The localized 
modes are centered around the three minima of a three-dimensional trapping potential for the BEC. The nonlinear 
interactions arise due to atomic collisions and conserve boson number. The coefficients for the nonlinear interactions 
are Xi> f° r sen? interactions, and X2, f° r intermodal collisions. These parameters depend on the s-wave scattering 
length for the condensate and its density as well as the degree of overlap between modes for xi- These coefficients of 
nonlinear dynamical terms are time-dependent: the scattering length can be varied by applying an external magnetic 
field to exploit a Feshbach resonance in the collisions between the atoms Jl7|. 

The time-varying linear quantum tunneling terms for the BEC, with coefficients Cljk, may be varied by modifying 
the optical potential barrier between minima of the trapping potential. The terms are left on for a short time to 
produce the desired SU(3) coherent state and then shut off. Of course they are not completely eliminated as some 
overlap is needed for the cross-mode collisions (with coefficient X2) to take place. It is important that the collision 
terms be strong enough to be nonnegligible even when linear quantum tunneling may be ignored. Tuning the scattering 
length is advantageous to ensure that the intermodal collision term is sufficiently strong. Although X2 might be weak, 
our subsequent analysis is valid for any choice of Xi ■ 

A typical initial condition could be N bosons in one region or mode. Allowing the system to evolve under the 
Hamiltonian ([!]) with Xi — 0, causes the system to evolve linearly, via coherent quantum tunneling, to a Perelomov 
SU(3) coherent state|| . As a special case, relevant to weak force detection |L3j and particle interferometry 113 , extremal 
coherent states are states which are mapped to zero under a ladder operator. For example, the SU(2) coherent 
state \j j) is annihilated by the specified raising operator J + , and the state \j — j) is annihilated by the specified 
lowering operator </_. These two states are extremal, and a superposition of these two extremal states would be 
(\j j) + \j Extremal SU(3) coherent states are states which can be annihilated by the su(3) ladder operators. 

From conservation of the total number N of bosons, the Hamiltonian may be rewritten in terms of the SU(3) 
generators and the total number N is the irreducible representation (irrep) parameter. The eight generators of SU(3) 
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may be defined as the two Cartan operators 

Xi = c\ci- 4c 2) (2) 



and the six generators 



X 2 = -(c\c 1+ clc 2 ~2clcl) (3) 

Y k = iiclcj - ctcfe) (4) 

Z k = clcj+c]c k , (5) 
where k = 1,2,3 and j = k mod 3 + 1. The raising operators and lowering operators are defined as 

4 = c{c 3 , (k<j), (6) 



Jl = (k>j), (7) 



or, alternatively, for k < j 



J% = \(Z k -iY k ) (8) 
Jl = \{Z k + iY k ), (9) 

in terms of SU(3) generators. 

The term involving the irrep parameter N can be removed by moving to a rotating picture. Furthermore, by 
assuming that the initial SU(3) coherent state has been generated from the evolution ([j]), with large linear quan- 
tum tunneling and negligible contributions from nonlinear evolution, the optical potential barrier between trapping 
potential minima is increased to reduce the linear tunneling. At the same time, the nonlinear coefficients must be 
reasonably large compared to the coefficient for linear tunneling, perhaps by employing an external magnetic field to 
exploit a Feshbach resonance, as discussed above. For large nonlincarities and small linear quantum tunneling terms, 
the Hamiltonian ([j]) can be approximated by 

H=^(X? + 3Xl), (10) 

with x = Xi ~Xi- The Hamiltonian is a sum of quadratic forms of Cartan operators which commute, and this property 
will be useful in later calculations. 



COHERENT STATES 
SU(2) coherent states 

In order to introduce the SU(3) coherent states, it is useful to review the SU(2) coherent states and to use this 
knowledge to generalize to SU(3) coherent states. SU(2) coherent states and their methods were first developed as 
atomic coherent states || [lO] to treat atomic systems. On the other hand, coherent states for Lie groups were defined 
by PerelomovQ as an orbit generated by group action on a reference state, which is usually the highest- (or lowest-) 
weight state. For compact groups, the simplest nontrivial case is SU(2), which can be parametrized by three real 
parameters {9, ipi, tpi}. The lowest order faithful representation of an arbitrary g 6 SU(2) is as a 2 x 2 matrix 

9(<PiA<f») = [_ e -^^e encase)' 

As SU(2) C SU(3), the three generators of SU(2) can be obtained from (g), (g), @ and (g). For example, the 
subgroup SU(2)i2 is generated by {Xi,Yi,Z{\. The Casimir invariant is J 2 = X\ + Yj 2 + Z\ with eigenvalue j(j + 1). 
The 2x2 representation corresponds to j = 1/2. The weight basis corresponds to \jm) with J z \jm) = m\jm). The 
highest -weight state is \jj) for which J+\jj) = 0. 

The SU(2) coherent state, for fixed j, is given by 



\9, tp) = exp 



lii). (12) 



The coherent state can be represented geometrically on the Poincare sphere, which is isomorphic to the coset space 
SU(2)/U(1). 
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SU(3) coherent states 

The SU(3) coherent states may be formulated as a generalization of the SU(2) coherent statesjl^, [lij, and are 
obtained by the action of SU(3) on the highest-weight state. In order to obtain the SU(3) coherent states, it is 
convenient to employ the decomposition of SU(3) in Ref. ^0|, whereby the SU(3) operator is decomposed into a 
combination of SU(2) operators. The decomposition allows us to parameterize an arbitrary element g in the 3x3 
matrix representation as 

/ 1 \ / e^costf -sin^ \ / 1 \ 

(13) 

where V(ipi, £, ip 2 ) and W(tf3, C, <fi) are 2x2 matrices representing elements of SU(2)23- The middle matrix on the 
right-hand side of Eq. ( |l3| ) is an element of SU(2)i2. 

The SU(3) generators of Eqs. (Q,||) and (^,||) are presented as a three-boson realization. It is convenient to employ 
the basis states, which are eigenstatcs of the SU(3) Casimir operators and of the two elements of the Cartan subalgebra 
(||,||). The basis state can be labeled by the three numbers n\, n 2 and n 3 , which satisfy 

n k \ni,n 2 ,n 3 ) = n k \ni, n 2 , n 3 ), for k € {1,2,3} . (14) 

As any action of the raising operators on the state \N, 0, 0) annihilates it, the state \N, 0, 0) is a highest-weight state. 
Coherent states are obtained by SU(3) action on this highest-weight state. 

Taking the highest-weight state \N, 0, 0) as the reference state, the action of Eq. ( |l3| ) on the highest-weight state 
gives an explicit form of coherent states. The factorization ( |l3| ) ensures that the matrix on the furthest right leaves the 
highest-weight state invariant, and only the two other matrices are important in determining the coherent state. The 
highest-weight state is thus invariant under SU(2)23; hence the coherent states are parametrized on the coset space 
SU(3)/SU(2). In fact the SU(3) coherent states can be parametrized on the coset space SU(3)/U(2) by eliminating 
an arbitrary phase, but the following calculations are easier if we retain the phase. In calculations of the Q~symbol 
for the state, this phase is not important. 

The coherent states can be generated as 

|£, 6, ip, ipi,(p 2 ) = g(<pi,£, <P2, <P, 0, <P3,(, <Pi)W, o, 0), (15) 

which can be expressed as 

/ AT\ 1/2 



31=0 V / 

3i , v 1/2 

x J2 e^V^'^W 2 ^) cos (jl - h \0 \N-jiJi (16) 

h=o ^ 32 ' 



This parameterization of the SU(3) coherent state includes an arbitrary phase. 



QUANTUM DYNAMICS OF THE SYSTEM 
Analytical solutions 

Let the initial state |^(0)) be an arbitrary SU(3) coherent state. The state \ip(t)} for arbitrary time t may thus be 
given as 

N / N \ 1/2 H f n \ 1/2 

|V(£)> = e MN - jl) sin jl (e) cos iN - jl) (6) J] e ih<e*j»pi(ji-i*) s in J ' 2 (f) cos^-^^f) Jl 

3i=o V 31 ' 32=0 V 32 ' 

xexp[-2ixt(N 2 /3 + j 2 +j 2 - 3l (N + j 2 ))] \N - j u j x - j 2 , j 2 ) , (17) 

under the Hamiltonian evolution (|l0|). The time-dependent element of Eq. (|l7]) exhibits the periodic evolution of the 
state and suggests that the state evolves into superposition states periodically, by analogy with the SU(2) coherent 
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state superposition caseQ, but with additional complexity due to the greater number of parameters and the higher 
dimension of the group manifold. 

From Eq. (|lj), we observe that the recurrence time, for which the state evolves cyclically back into the original 
state, is given by r = 7r/x for all values of the irrep parameter N. At half the recurrence time, r/2 = tt/2x, the state 
evolves into the superposition 

|V(r/2)) - i e -^ 2 /3 07 ^ VuV2) + ^ 0, ^ Vl + 7F] V2) + ^ ^ ^ ^ ; + ^ ^ ^ Vl + 7r ^ 2 + n )] . (is) 
for even N , or 

|V(r/2)) = l e -^ 2 /3 flj ^ ^ ^ + |^ ^ + 7rj ^^ 2 ) + |£ ^ ^ ^ + 7r? + |^ 0; ^ + „■)] , (19) 

for odd N. At one-quarter of the recurrence time r/4, the state evolves into a superposition state with the exact 
form depending on the total number N. There are four types of superposition states at r/4, classified by the quantity 
N mod 4. The four types are given as the following, 



1. If N — An, for n an integer, 
1 



rl> (r/4)) = -e 



-inN 2 /6 



|£, 9, Lp, pi, p 2 ) + |£, 9, ip, (pi + n, <p 2 ) + |£, 9, ip, tpi,ip2+n) + |£, 0, tp, ipx +ir,tp 2 + it) 



+ i 



7T \ 7T 
£, 6, Ifi, (ft, lfi 2 - 2/ _ l& 0, Pi . Pli P2 + - ) 



7T 7T 

£,0,P,Pi + 2'^ 2 - 2 



7T 7T 

£,6,tp,ipl + -,ip 2 + 2 



7T 7T 7T 7T 

|£,0,p,Pl + -,<^2 +7T) - |£,0,P,Pl + - ,Pa) + |£,0,p,Pl +TT,<p 2 - ^) + I^J ^) Vl +n,lp2 + ~) 

TT TT 7T 7T TT 7T 

|£, 19, p, ipi - -, lp 2 - -) + |£, 9, (f, tpi - -, (f 2 + -) - |£, 9, Lp, Ifl - -, Lp 2 ) + |£, 9,Lp,L Pl - -,Lp 2 + 7T) 



20) 



2. If N = 4n + 1, 



-i-KN 2 /6 



|£, 0, V) Pl + 2 , P2 + ^ ) + |£, 0, P, Pi ~ ~, P3 + ~) + l£) P) Pl + 2 > P2 - ^> + l£> 0) Pi Pi ~ 2 > P2 - 2 ) 

/ 7T 7T 

+ M - If, 0) Pi Pi + 2 ' ^ 2 > - l£) 0, Pj Pi + 2 ' ^2 + ^ + l£' 0' V' ^1 + 7F ' ^2) ~ |£, 0, P, Pi + 7T, P2 + 7T> 

7T 7T 7T 7T 

+ |£) 0, P, Pl + 1", P2 " 77> - Ifi 0, Pi Pi + T) p 2 + 77) + l£i ^i Pi Pi - o 1 P2> + ICi 0i Pi Pi - 77i P2 + T> 



7T 7T 
ICi 0i Pi Pli P2) + ICi 0i Pi PlP2 + 7T) - |C, 0, P, Pl, P2 + -g) + |£) 0, P, Pl, Pa - g ) 



(21) 



3. If TV = 4n + 2, 



^ (r/4)) = -e 



-iirN 2 /6 



|£i 0i Pi Pi, P2) + |Ci 0, Pi Pi + 7T) Pa) + ICi 0i Pi Pl, Pa + tt) + |^, 9, ip, ipx +ir,Lp 2 + tt) 

(7T TT TT 7T TT TT 

- ICi 0, Pi Pl + ""i P2 + ^ ) _ 16 0' P' Pl + ^i P2 - 2) + 1^' ' V' ^1 ~ 2 ' ^ 2 + 2^ _ </ '' V ' 1 ~ 2 ,ip2 ~ 2 ^ 

7T 7T 7T 7T 

+ ICi0iPiPl - ~ 1^'0'P'Pl - 2'^ 2 + l£,0iP,Pl,P2 + g ) + l£>0iPiPliP2 - 2^ 



TT TT TT TT TT TT 

|£, 0, P, Pi + 2 ' ^2 + 2 ) + 0' P' Pi + 2 ' P2 - 2^ ~ <y9 ' + 2 ,<P2+n ^ + ^' 6 ' if ' ¥?1 + 2 ' 



(22) 
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4. If TV = An + 3, 



I 1> (r/4)) = i e -^V6 



7T7T 7T7T 7T7T 7T7T 

|£, V) <Pl ~ 2 > </>2 - - ) + |£, 0, <f, Ifl + -,ip 2 - -) + |£, 6, if, (fix - -, f 2 + -> + |£, 0, ¥>, + 2 ' ^ 2 + 2 ' 

/ 7T 7T 

7T 7T 7T 7T 



7T 7T 
|£, 0, i/3, if i + IT, If 2 + 7T) + |£, 0, V, + 7T, 99 2 ) - |£, 0, ¥>1 + ""j ¥>2 - g ) + ^ ¥>> +Tr,<p 2 + -) 



(23) 



The Q— symbol 

In the previous section we have obtained analytical expressions of \^>{t)) for two particular time values, r/2 and 
r/4. More generally however we can obtain the state analytically for all values of t. The density matrix pit") for the 
system may be written in terms of the pure state \i/)(t)) as 



The Q-symbol for this state is given 



P (t) = m))(m\- 



Q(t) = (£, o, if, <pi,v>2\p(t)\€, o, f, if 1 ,if 2 )- 



(24) 



(25) 



The arbitrary phase disappears in the calculation of Q(t). It is thus convenient to employ the two relative phases 
4>i = fx ~~ f an d 4>2 = f2 — fi ■ The Q-symbol for the state may thus be expressed as 



Q^,e,<t>i,<h-m,m,M0),M0);t) 



N h 
01 =0j 2 =0 



(26) 



where 



h ,32 



-in -0! (o)) / gin g gin g ^y 1 (- cos g CQS 9 ( Q j 



,N- n ( N 



Jl 

itf-jl / Jl 



x e -w^-^( ))(sinCsm^(0)) il (cos^cos^(0))' - j v ] < 



(27) 



The Q-symbol depends on four periodic parameters, and therefore a plot of the Q-symbol would need to be 
embedded in five-dimensional space to be viewed. However, in order to view the plot we need to reduce the number 
of dimensions required for plotting the function. A simple way to plot the Q-symbol is to fix one parameter, for 
example by setting 6 — ir/2. This corresponds to a slice of the Q-symbol, and the choice of where to slice depends 
on the features to be emphasized in the plot. A slice is a standard means for representing multi-dimensional figures; 
animation can be used, for example, to present the figure by 'traveling' through a sequence of slices. We will present 
an example in the following subsection. 



An Example 

The nonlinear SU(2) oscillator generates superposition states from an initial SU(2) coherent state and can generate 
a superposition of two distinct SU(2) coherent states at a half of its recurrence time ||. This superposition of two 
distinct SU(2) coherent states results in two major peaks arising in the plot of the Q-symbol for the superposition 
state. Interference fringes midway between the two peaks arise due to the coherence properties of the superposition. 
These interference fringes can be clearly seen as dimples in a plot of the logarithm of the Q-symbol]2l|. 

The superposition states we have seen in Eqs. (|20|-p3|) arise from the nonlinearity of the SU(3) system. In this 
subsection, we wish to observe a manifestation of such a superposition by plotting the Q-symbol for the state. As 
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discussed in the previous subsection, we can reduce the dimension of the parameter space to two by obtaining a slice 
of the Q-symbol for the state. We consider an SU(3) coherent state which is also an SU(2) 2 3 coherent state. This 
restriction to SU(2) 2 3 coherent states ensures that only a subset of the Q-symbol is explored, with this subset being a 
slice which can be plotted on the Poincare sphere corresponding to the SU(2) 2 3/U(1) coset space. An initial SU(2) 2 3 
coherent state can be written as 

|^(0)) = |£(0),tt/2,<M0),<M0)) 

N / N \ 1/2 

= e -#i(o)JV y- e iMa)h sin h^f \\w S ( N -i 2 \£(0)) [. \0,N-j 2 ,j 2 ) (28) 

A state in the SU(2) 2 3 subsystem stays in this subsystem under time evolution, so the factor of <pi becomes an 
arbitrary phase and disappears in the calculation of the Q -s ymbol. We now employ the notation Q to refer to a slice 
of the Q-symbol. In this notation, the reduced Q-symbol, with the initial condition (|28|), is given by 

Q(0,£,0 2 ;i) = ^,e,<h,Mp(fMM,<h) 

N N 

\P+<7/ ,t\ „ /^n\\\( 2Ar -(P+9)) 



sm 2N (6) e-***-^-*'*™ ( sin(e) sin(£(0))) P+<? ( cos(£) cos(£(0))) ( 

:)(: 

2N, 



X 



exp[-2i X t(p-q)( P -q-N)} . (29) 



The ^-component of the Q -s ymbol i§ time-independent with a fixed weighting of sin (9). The optimal slice for 
observing the Q-symbol in the reduced £-^ 2 space is thus obtained by setting 9 = ir/2. 
With the initial condition £(0) = 7r/4, </> 2 (0) = 0, Q is simplified to 



N N 



Q(£,0 2 : t) = 2- N Y,H z- i(p - 9) * 2 ( " ) ( N a ) sin^(£) cos^"^ (£) 



p=Q q=0 



P J \ 1 

xexp[-2i X t(p-q)(p + q-N)} , (30) 



and is numerically represented in Fig. [l] and Fig. ^| The time-dependent term in Eq. ( p0| ) implies that the time- 
evolution of the state is distinguishable in terms of the total number N. The product (p — q) (p + q — N) is even for any 
p and q if N is odd, whereas it is not true for any even total number N. The recurrence time of the time evolution 
with an odd total number N is a half of the recurrence time with an even total number. We solve this evolution 
numerically and plot Q as it evolves for two values of the total number, N = 9, 10. Fig. |l| shows the time-evolution 
of Q for the state for even number N = 10; in contrast, Fig. § depicts the time-evolution of the Q for the state for 
odd number N — 9. 



CONCLUSIONS 



We have investigated the formation of superpositions of SU(3) coherent states via a nonlinear Hamiltonian evolution. 
This evolution could be realized in terms of four-wave mixing in nonlinear optics or in terms of a nonlinear interaction 
between three independent modes of a Bose-Einstein condensate with nonncgligiblc nonlinear interactions between 
the separate, but overlapping, modes. The linear quantum tunneling term is allowed to vary independently of the 
nonlinear (collision) terms, and linear quantum tunneling is used to prepare the SU(3) coherent state for iV bosons 
in the three modes. Then this SU(3) coherent state undergoes a nonlinear evolution into a superposition of SU(3) 
coherent states. 

We have obtained explicit expressions for the superposition of coherent states at one-half and one-quarter the 
recurrence time for this periodic evolution. Whereas superpositions of pairs of coherent states are obtained for half 
the recurrence time for Heisenberg-Weyl, SU(2) and SU(1,1) systems, the higher dimension of SU(3) dynamics leads to 
much more complex and interesting superpositions even at half the recurrence time. The nature of the superposition 
depends on the quantity N mod 4, with N being the total quantum number. The state at some time t can be 
represented by the use of a Q-symbol representation. However, the multidimensional domain of the Q-symbol makes 
the visualization of this function challenging. We have suggested that the method of plotting slices of the Q-symbol is 
desirable in this respect. An example was provided where the dynamics could be restricted to SU(2) 2 3, which makes 
the plots quite clear and readily interpreted: a superposition of two SU(2) coherent states is evident in these plots. 
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t=0 t=0. 19635 t=0.3926S9 




t=0 589049 t=0. 785398 t=0.9B1748 




t-1.1781 t-1.37445 t-1.5708 




FIG. 1: Q(€,,fa) is plotted at a certain time t with the use of the stereographical mapping. The origin corresponds the north 
pole of the spherical subdomain (£, fa). Time flows from the left top figure at t = to the right bottom figure at t — r/2. The 
recurrence time is twice as longer as the odd number case below. 



t=0 t=0. 19635 t=0.392699 




t-0.E8904S t-0.7853 98 t-0.S81748 




t=1 1781 t=1.37445 t=1.5708 




FIG. 2: Q(^,fa) is plotted at a certain time t with the use of the stereographical mapping. The origin corresponds the north 
pole of the spherical subdomain (£, fa). Time flows from the left top figure at t = to the right bottom figure at t — r/2. The 
recurrence time is a half of the even number case above. 

This approach may be extended to SU(n) coherent states by generalizing the Hamiltonian (Q) to n interacting 
modes. However, the SU(3) dynamics is the most relevant case to current physical realizations. Although nonlinear 
optics and Bose-Einstein condensation have been specifically mentioned, any three-boson system with the given 
nonlinear evolution (|l]) could yield superpositions of SU(3) coherent states, provided that the initial state is itself an 
SU(3) coherent state. 
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